Identification of Potential Long Non-Coding RNA Candidates that Contribute to Triple-Negative Breast Cancer in Humans through Computational Approach.

Breast cancer (BC) is the most frequent malignancy identified in adult females, resulting in enormous financial losses worldwide. Owing to the heterogeneity as well as various molecular subtypes, the molecular pathways underlying carcinogenesis in various forms of BC are distinct. Therefore, the advancement of alternative therapy is required to combat the ailment. Recent analyses propose that long non-coding RNAs (lncRNAs) perform an essential function in controlling immune response, and therefore, may provide essential information about the disorder. However, their function in patients with triple-negative BC (TNBC) has not been explored in detail. Here, we analyzed the changes in the genomic expression of messenger RNA (mRNA) and lncRNA in standard control in response to cancer metastasis using publicly available single-cell RNA-Seq data. We identified a total of 197 potentially novel lncRNAs in TNBC patients of which 86 were differentially upregulated and 111 were differentially downregulated. In addition, among the 909 candidate lncRNA transcripts, 19 were significantly differentially expressed (DE) of which three were upregulated and 16 were downregulated. On the other hand, 1901 mRNA transcripts were significantly DE of which 1110 were upregulated and 791 were downregulated by TNBCs subtypes. The Gene Ontology (GO) analyses showed that some of the host genes were enriched in various biological, molecular, and cellular functions. The Kyoto encyclopedia of genes and genomes (KEGG) pathway analysis showed that some of the genes were involved in only one pathway of prostate cancer. The lncRNA-miRNA-gene network analysis showed that the lncRNAs TCONS_00076394 and TCONS_00051377 interacted with breast cancer-related micro RNAs (miRNAs) and the host genes of these lncRNAs were also functionally related to breast cancer. Thus, this study provides novel lncRNAs as potential biomarkers for the therapeutic intervention of this cancer subtype.

RMVar: an updated database of functional variants involved in RNA modifications.

Distinguishing the few disease-related variants from a massive number of passenger variants is a major challenge. Variants affecting RNA modifications that play critical roles in many aspects of RNA metabolism have recently been linked to many human diseases, such as cancers. Evaluating the effect of genetic variants on RNA modifications will provide a new perspective for understanding the pathogenic mechanism of human diseases. Previously, we developed a database called 'm6AVar' to host variants associated with m6A, one of the most prevalent RNA modifications in eukaryotes. To host all RNA modification (RM)-associated variants, here we present an updated version of m6AVar renamed RMVar (http://rmvar.renlab.org). In this update, RMVar contains 1 678 126 RM-associated variants for 9 kinds of RNA modifications, namely m6A, m6Am, m1A, pseudouridine, m5C, m5U, 2'-O-Me, A-to-I and m7G, at three confidence levels. Moreover, RBP binding regions, miRNA targets, splicing events and circRNAs were integrated to assist investigations of the effects of RM-associated variants on posttranscriptional regulation. In addition, disease-related information was integrated from ClinVar and other genome-wide association studies (GWAS) to investigate the relationship between RM-associated variants and diseases. We expect that RMVar may boost further functional studies on genetic variants affecting RNA modifications.

Lnc2Cancer 3.0: an updated resource for experimentally supported lncRNA/circRNA cancer associations and web tools based on RNA-seq and scRNA-seq data.

An updated Lnc2Cancer 3.0 (http://www.bio-bigdata.net/lnc2cancer or http://bio-bigdata.hrbmu.edu.cn/lnc2cancer) database, which includes comprehensive data on experimentally supported long non-coding RNAs (lncRNAs) and circular RNAs (circRNAs) associated with human cancers. In addition, web tools for analyzing lncRNA expression by high-throughput RNA sequencing (RNA-seq) and single-cell RNA-seq (scRNA-seq) are described. Lnc2Cancer 3.0 was updated with several new features, including (i) Increased cancer-associated lncRNA entries over the previous version. The current release includes 9254 lncRNA-cancer associations, with 2659 lncRNAs and 216 cancer subtypes. (ii) Newly adding 1049 experimentally supported circRNA-cancer associations, with 743 circRNAs and 70 cancer subtypes. (iii) Experimentally supported regulatory mechanisms of cancer-related lncRNAs and circRNAs, involving microRNAs, transcription factors (TF), genetic variants, methylation and enhancers were included. (iv) Appending experimentally supported biological functions of cancer-related lncRNAs and circRNAs including cell growth, apoptosis, autophagy, epithelial mesenchymal transformation (EMT), immunity and coding ability. (v) Experimentally supported clinical relevance of cancer-related lncRNAs and circRNAs in metastasis, recurrence, circulation, drug resistance, and prognosis was included. Additionally, two flexible online tools, including RNA-seq and scRNA-seq web tools, were developed to enable fast and customizable analysis and visualization of lncRNAs in cancers. Lnc2Cancer 3.0 is a valuable resource for elucidating the associations between lncRNA, circRNA and cancer.

LnCeCell: a comprehensive database of predicted lncRNA-associated ceRNA networks at single-cell resolution.

Within the tumour microenvironment, cells exhibit different behaviours driven by fine-tuning of gene regulation. Identification of cellular-specific gene regulatory networks will deepen the understanding of disease pathology at single-cell resolution and contribute to the development of precision medicine. Here, we describe a database, LnCeCell (http://www.bio-bigdata.net/LnCeCell/ or http://bio-bigdata.hrbmu.edu.cn/LnCeCell/), which aims to document cellular-specific long non-coding RNA (lncRNA)-associated competing endogenous RNA (ceRNA) networks for personalised characterisation of diseases based on the 'One Cell, One World' theory. LnCeCell is curated with cellular-specific ceRNA regulations from >94 000 cells across 25 types of cancers and provides >9000 experimentally supported lncRNA biomarkers, associated with tumour metastasis, recurrence, prognosis, circulation, drug resistance, etc. For each cell, LnCeCell illustrates a global map of ceRNA sub-cellular locations, which have been manually curated from the literature and related data sources, and portrays a functional state atlas for a single cancer cell. LnCeCell also provides several flexible tools to infer ceRNA functions based on a specific cellular background. LnCeCell serves as an important resource for investigating the gene regulatory networks within a single cell and can help researchers understand the regulatory mechanisms underlying complex microbial ecosystems and individual phenotypes.

RMDisease: a database of genetic variants that affect RNA modifications, with implications for epitranscriptome pathogenesis.

Deciphering the biological impacts of millions of single nucleotide variants remains a major challenge. Recent studies suggest that RNA modifications play versatile roles in essential biological mechanisms, and are closely related to the progression of various diseases including multiple cancers. To comprehensively unveil the association between disease-associated variants and their epitranscriptome disturbance, we built RMDisease, a database of genetic variants that can affect RNA modifications. By integrating the prediction results of 18 different RNA modification prediction tools and also 303,426 experimentally-validated RNA modification sites, RMDisease identified a total of 202,307 human SNPs that may affect (add or remove) sites of eight types of RNA modifications (m6A, m5C, m1A, m5U, Ψ, m6Am, m7G and Nm). These include 4,289 disease-associated variants that may imply disease pathogenesis functioning at the epitranscriptome layer. These SNPs were further annotated with essential information such as post-transcriptional regulations (sites for miRNA binding, interaction with RNA-binding proteins and alternative splicing) revealing putative regulatory circuits. A convenient graphical user interface was constructed to support the query, exploration and download of the relevant information. RMDisease should make a useful resource for studying the epitranscriptome impact of genetic variants via multiple RNA modifications with emphasis on their potential disease relevance. RMDisease is freely accessible at: www.xjtlu.edu.cn/biologicalsciences/rmd.

Construction and investigation of a combined hypoxia and stemness index lncRNA-associated ceRNA regulatory network in lung adenocarcinoma.

Hypoxia and stemness are important factors in tumor progression. We aimed to explore the ncRNA classifier associated with hypoxia and stemness in lung adenocarcinoma (LUAD). We found that the prognosis of LUAD patients with high hypoxia and stemness index was worse than that of patients with low hypoxia and stemness index. RNA expression profiles of these two clusters were analyzed, and 6867 differentially expressed (DE) mRNAs were screened. Functional analysis showed that DE mRNAs were associated with cell cycle and DNA replication. Protein-protein interaction network analysis revealed 20 hub genes, among which CENPF, BUB1, BUB1B, KIF23 and TTK had significant influence on prognosis. In addition, 807 DE lncRNAs and 243 DE miRNAs were identified. CeRNA network analysis indicated that AC079160.1-miR-539-5p-CENPF may be an important regulatory axis that potentially regulates the progression of LUAD. The expression of AC079160.1 and CENPF were positively correlated with hypoxia and stemness index, while miR-539-5p expression level was negatively correlated with hypoxia and stemness index. Overall, we identified CENPF, BUB1, BUB1B, KIF23 and TTK as potentially key genes involved in regulating hypoxia-induced tumor cell stemness, and found that AC079160.1-miR-539-5p-CENPF axis may be involved in regulating hypoxia induced tumor cell stemness in LUAD.

Predicting miRNA-disease association from heterogeneous information network with GraRep embedding model.

In recent years, accumulating evidences have shown that microRNA (miRNA) plays an important role in the exploration and treatment of diseases, so detection of the associations between miRNA and disease has been drawn more and more attentions. However, traditional experimental methods have the limitations of high cost and time- consuming, a computational method can help us more systematically and effectively predict the potential miRNA-disease associations. In this work, we proposed a novel network embedding-based heterogeneous information integration method to predict miRNA-disease associations. More specifically, a heterogeneous information network is constructed by combining the known associations among lncRNA, drug, protein, disease, and miRNA. After that, the network embedding method Learning Graph Representations with Global Structural Information (GraRep) is employed to learn embeddings of nodes in heterogeneous information network. In this way, the embedding representations of miRNA and disease are integrated with the attribute information of miRNA and disease (e.g. miRNA sequence information and disease semantic similarity) to represent miRNA-disease association pairs. Finally, the Random Forest (RF) classifier is used for predicting potential miRNA-disease associations. Under the 5-fold cross validation, our method obtained 85.11% prediction accuracy with 80.41% sensitivity at the AUC of 91.25%. In addition, in case studies of three major Human diseases, 45 (Colon Neoplasms), 42 (Breast Neoplasms) and 44 (Esophageal Neoplasms) of top-50 predicted miRNAs are respectively verified by other miRNA-disease association databases. In conclusion, the experimental results suggest that our method can be a powerful and useful tool for predicting potential miRNA-disease associations.

Dysregulation of miRNAs as a signature for diagnosis and prognosis of gastric cancer and their involvement in the mechanism underlying gastric carcinogenesis and progression.

Gastric cancer (GC) is the second main cause of cancer-related mortality worldwide. The poor prognosis and survival of GC are due to diagnosis in an advanced, noncurable stage and with a restricted response to chemotherapy. GC is usually monitored in an advanced stage; therefore, the poor prognosis and lower level of survival rate with a restricted response to chemotherapy can be detected. Valuable and sensible biomarkers are urgently needed to display screen patients with a high risk of GC that can complement endoscopic diagnosis. Such biomarkers will enable the efficient prediction of the therapeutic response and prognosis of GC patients and prefer the establishment of an advantageous treatment method for each and every patient. Noninvasive diagnostic biomarkers may additionally make a contribution to the early identification of GC and enhance medical management. MicroRNAs (miRNAs) are a group of small noncoding RNAs that have displayed a strong association with GC. Accumulating evidence indicates that miRNAs are potential biomarkers with more than one diagnostic function for GC. Actually, miRNAs regulate cell proliferation, apoptosis, migration, invasion, and metastasis via many biological pathways through the repression of target mRNAs. The current review is accordingly to spotlight the multifaceted roles of miRNAs in GC, which would provide indications for future research. Therefore, we review right here the aberrant expression of miRNAs and underlying mechanisms, consequent effects due to miRNAs dysregulation, and accountable target genes in GC. Besides, potential clinical applications are also highlighted.

microRNA classifiers are powerful diagnostic/prognostic tools in ALK-, EGFR-, and KRAS-driven lung cancers.

microRNAs (miRNAs) can act as oncosuppressors or oncogenes, induce chemoresistance or chemosensitivity, and are major posttranscriptional gene regulators. Anaplastic lymphoma kinase (ALK), EGF receptor (EGFR), and V-Ki-ras2 Kirsten rat sarcoma viral oncogene homolog (KRAS) are major drivers of non-small cell lung cancer (NSCLC). The aim of this study was to assess the miRNA profiles of NSCLCs driven by translocated ALK, mutant EGFR, or mutant KRAS to find driver-specific diagnostic and prognostic miRNA signatures. A total of 85 formalin-fixed, paraffin-embedded samples were considered: 67 primary NSCLCs and 18 matched normal lung tissues. Of the 67 primary NSCLCs, 17 were echinoderm microtubule-associated protein-like 4-ALK translocated (ALK(+)) lung cancers; the remaining 50 were not (ALK(-)). Of the 50 ALK(-) primary NSCLCs, 24 were EGFR and KRAS mutation-negative (i.e., WT; triple negative); 11 were mutant EGFR (EGFR(+)), and 15 were mutant KRAS (KRAS(+)). We developed a diagnostic classifier that shows how miR-1253, miR-504, and miR-26a-5p expression levels can classify NSCLCs as ALK-translocated, mutant EGFR, or mutant KRAS versus mutation-free. We also generated a prognostic classifier based on miR-769-5p and Let-7d-5p expression levels that can predict overall survival. This classifier showed better performance than the commonly used classifiers based on mutational status. Although it has several limitations, this study shows that miRNA signatures and classifiers have great potential as powerful, cost-effective next-generation tools to improve and complement current genetic tests. Further studies of these miRNAs can help define their roles in NSCLC biology and in identifying best-performing chemotherapy regimens.

A Systematic Evaluation of Feature Selection and Classification Algorithms Using Simulated and Real miRNA Sequencing Data.

Sequencing is widely used to discover associations between microRNAs (miRNAs) and diseases. However, the negative binomial distribution (NB) and high dimensionality of data obtained using sequencing can lead to low-power results and low reproducibility. Several statistical learning algorithms have been proposed to address sequencing data, and although evaluation of these methods is essential, such studies are relatively rare. The performance of seven feature selection (FS) algorithms, including baySeq, DESeq, edgeR, the rank sum test, lasso, particle swarm optimistic decision tree, and random forest (RF), was compared by simulation under different conditions based on the difference of the mean, the dispersion parameter of the NB, and the signal to noise ratio. Real data were used to evaluate the performance of RF, logistic regression, and support vector machine. Based on the simulation and real data, we discuss the behaviour of the FS and classification algorithms. The Apriori algorithm identified frequent item sets (mir-133a, mir-133b, mir-183, mir-937, and mir-96) from among the deregulated miRNAs of six datasets from The Cancer Genomics Atlas. Taking these findings altogether and considering computational memory requirements, we propose a strategy that combines edgeR and DESeq for large sample sizes.

Discovery of novel splice forms and functional analysis of cancer-specific alternative splicing in human expressed sequences.

We report here a genome-wide analysis of alternative splicing in 2 million human expressed sequence tags (ESTs), to identify splice forms that are up-regulated in tumors relative to normal tissues. We found strong evidence (P < 0.01) of cancer-specific splice variants in 316 human genes. In total, 78% of the cancer-specific splice forms we detected are confirmed by human-curated mRNA sequences, indicating that our results are not due to random mis-splicing in tumors; 73% of the genes showed the same cancer-specific splicing changes in tissue-matched tumor versus normal datasets, indicating that the vast majority of these changes are associated with tumorigenesis, not tissue specificity. We have confirmed our EST results in an independent set of experimental data provided by human-curated mRNAs (P-value 10(-5.7)). Moreover, the majority of the genes we detected have functions associated with cancer (P-value 0.0007), suggesting that their altered splicing may play a functional role in cancer. Analysis of the types of cancer-specific splicing shifts suggests that many of these shifts act by disrupting a tumor suppressor function. Sur prisingly, our data show that for a large number (190 in this study) of cancer-associated genes cloned originally from tumors, there exists a previously uncharacterized splice form of the gene that appears to be predominant in normal tissue.

Decay rates of human mRNAs: correlation with functional characteristics and sequence attributes.

Although mRNA decay rates are a key determinant of the steady-state concentration for any given mRNA species, relatively little is known, on a population level, about what factors influence turnover rates and how these rates are integrated into cellular decisions. We decided to measure mRNA decay rates in two human cell lines with high-density oligonucleotide arrays that enable the measurement of decay rates simultaneously for thousands of mRNA species. Using existing annotation and the Gene Ontology hierarchy of biological processes, we assign mRNAs to functional classes at various levels of resolution and compare the decay rate statistics between these classes. The results show statistically significant organizational principles in the variation of decay rates among functional classes. In particular, transcription factor mRNAs have increased average decay rates compared with other transcripts and are enriched in "fast-decaying" mRNAs with half-lives <2 h. In contrast, we find that mRNAs for biosynthetic proteins have decreased average decay rates and are deficient in fast-decaying mRNAs. Our analysis of data from a previously published study of Saccharomyces cerevisiae mRNA decay shows the same functional organization of decay rates, implying that it is a general organizational scheme for eukaryotes. Additionally, we investigated the dependence of decay rates on sequence composition, that is, the presence or absence of short mRNA motifs in various regions of the mRNA transcript. Our analysis recovers the positive correlation of mRNA decay with known AU-rich mRNA motifs, but we also uncover further short mRNA motifs that show statistically significant correlation with decay. However, we also note that none of these motifs are strong predictors of mRNA decay rate, indicating that the regulation of mRNA decay is more complex and may involve the cooperative binding of several RNA-binding proteins at different sites.